Mult i— scalar black holes with contingent primary hair: Mechanics and stability 
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We generalize a class of magnetically charged black holes holes non-minimally coupled to two 
scalar fields previously found by one of us to the case of multiple scalar fields. The black holes 
possess a novel type of primary scalar hair, which we call a contingent primary hair, although the 
solutions possess degrees of freedom which are not completely determined by the other charges of the 
theory, the charges necessarily vanish in the absence of the magnetic monopole. Only one constraint 
relates the black hole mass to the magnetic charge and scalar charges of the theory. We obtain a 
Smarr-type thermodynamic relation, and the first law of black hole thermodynamics for the system. 
We further explicitly show in the two-scalar-field case that, contrary to the case of many other 
hairy black holes, the black hole solutions are stable to radial perturbations. 
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I. INTRODUCTION 

In classical general relativity, no hair theorems im- 
pose strong constraints on the possibility of obtaining 
black hole solutions of the Einstein equations coupled to 
non-trivial scalar fields. A crucial ingredient for their 
proof is that the scalars be minimally coupled to grav- 
ity and other fields. When this condition is relaxed new 
possibilities emerge for evading the no hair theorems. 

One of the first investigations was undertaken by 
Bekenstein [3] , who attempted to show that a non-trivial 
black hole solution exists for a conformally coupled scalar 
field. Although the scalar field diverges at the horizon, 
Bekenstein argued that this solution nevertheless admits 
a physical interpretation. Unfortunately, a more detailed 
analysis of the problem casts serious doubt on this 
We will therefore take the view that scalar fields should 
be regular at the horizon for genuine black hole solutions. 

Non-trivial scalar hair is possible, however, when black 
holes are coupled to scalar fields with non-linear self- 
interactions. Such solutions were first found in the case 
of gravity coupled to the Skyrme model |^J, and subse- 
quently for the Einstein- Yang-Mills-Higgs model , and 
other generalizations. The large literature on this topic 
has been recently reviewed in ref. j^. What appears to 
be characteristic of the class of black holes with non- 
linear self-interactions is that the scalar fields fall off very 
rapidly at spatial infinity - e.g., exponentially fast - and 
hence the asymptotic scalar charges vanish. The scalar 
hair is characterized instead by non-trivial charges on the 
horizon. 
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A third possibility is that of minimally coupled scalar 
fields with potentials, V{(^, which violate the dominant 
energy condition (DEC) 0,1111. In single scalar field 
models, the DEC is equivalent to V{(j)) > 0, and thus 
one must choose potentials for which V{(j)) < for some 
field values. Examples of asymptotically flat black hole 
solutions have been found for cases in which ¥{(/)) < 
everywhere "s^, and for cases in which V{(jji) possesses at 
least one global minimum at negative values Jt,tM- Both 
analytic 0,13 and numerical 0,13 examples are known. 

A final possibility is to consider black holes with scalar 
fields which are non-minimally coupled to gauge fields. 
Such models have been extensively investigated because 
they arise naturally in Kaluza-Klein theories and in the 
effective low-energy limit of string theory, where the dila- 
ton plays a non-trivial role. In all these models one can 
find black hole solutions with non-zero scalar charges at 
spatial infinity |lOl [T]| . An analogous phenomenon can 
be shown to take place also in the pure dilaton-gravity 
sector of effective string theories, when one takes into 
account the coupling of the dilaton to gravity via Gauss- 
Bonnet terms [ij. 

For this class of non-minimally coupled models with 
non-zero asymptotic scalar charges, however, the scalar 
charges in question are not independent parameters, but 
in all cases are a given function of the other asymptotic 
charges which characterize the solution, namely the ADM 
mass and the electric and magnetic charges etc. As a re- 
sult such scalar charges have been called secondary hair 
by the authors of Refs. [Hill, to distinguish them from 
the theoretical possibility of a primary hair, namely an 
asymptotic scalar charge which is completely indepen- 
dent of the other charges. 

In light of all the known solutions ^ one may therefore 
conjecture 14] that a weaker form of the no- hair prop- 
erty still holds, namely that for theories which satisfy 
the DEC black hole solutions can be classified by a small 
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number of parameters, which in the case of asymptotic 
charges include only conserved charges such as mass, an- 
gular momentum and gauge charges, but no asymptotic 
scalar charges. 

Even this no hair conjecture does not appear to be 
valid in general, however. In a recent paper by one of us 
the properties of magnetically charged black holes 
coupled to a dilaton, and an additional modulus field, 
according the the field equations generated by the 4- 
dimensional action |16| 
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were studied using techniques from the general theory of 
dynamical systems, which have been previously applied 
to static spherically symmetric solutions of gravity cou- 
jled to scalar fields in a number of contexts [T^l • In Ref. 



it was shown that the regular black hole solutions 
were parameterized by an additional degree of freedom 
in addition to the mass, M, and magnetic charge, Q, and 
it was conjectured that this degree of freedom could be 
considered to be a "primary scalar hair" . 

An important issue in this context is that of stability 
of the hairy solutions. The physical relevance of the so- 
lutions would in fact be spoiled by the presence of insta- 
bilities. Instabilities are known to be present in the case 
of Einstein- Yang-Mills models and the DEC- 

violating solutions for which the stability issue has been 
most thoroughly studied 0. Dilaton black holes with 
a single scalar do not appear to suffer from instability 
problems [23ll2^: however, they possess only secondary 
hair. 

It is the aim of the present paper to clarify and extend 
the results of Ref. jTgl by determining various relations 
satisfied by the charges, and discussing the stability of 
the solutions. We will find in particular, that although 
the scalar degree of freedom cannot be considered to be 
a primary hair in the strictest sense, since it must nec- 
essarily vanish if the electromagnetic field vanishes, the 
regular static spherically symmetric black hole solutions 
are not completely specified by their mass and magnetic 
charge. The solutions therefore are potentially of consid- 
erable physical interest as the only known static spheri- 
cally symmetric solutions with non-trivial scalar charges 
at spatial infinity which are not completely determined 
by the other asymptotic charges of the theory, and there- 
fore provide a counter-example to some forms of the no- 
hair conjecture. We will show that the solutions are sta- 
ble, adding significance to their interpretation. 

Intuitively, one might attempt to understand the phys- 
ical basis of the no-hair conjecture as arising from the fact 
there is no scalar equivalent of Gauss's law to give a con- 
served charge. The presence of the electromagnetic field 



supports the scalar charge in the case of models with a 
secondary scalar hair since the scalar coupling enters into 
the equivalent of Gauss's law. In the solutions described 
here the scalar charges are also supported by the pres- 
ence of the electromagnetic field. However, they are not 
entirely determined by it. 

It is quite consistent with past usage to adopt the ter- 
minology "primary hair" for the additional degrees of 
freedom which arise in the model However, the fact 
that the scalar charges are not entirely independent of the 
electromagnetic field suggests a new terminology might 
be appropriate, as a way of capturing the finer distinc- 
tions that the present model has revealed. We will there- 
fore refer to the scalar charges which could theoretically 
exist in the absence of additional non-zero gauge charges 
as elementary primary hair, as would be consistent with 
the original aims of the first "no hair" theorems Q . The 
additional degrees of freedom which arise in the present 
model could by contrast be deemed to be a contingent 
primary hair. 



II. MULTI-SCALAR BLACK HOLE SOLUTIONS 

Rather than restricting our attention to solutions of the 
field equations obtained from varying the action JQ), we 
will instead consider the somewhat more general action 
for N scalar fields coupled to a single [/(I) Abelian gauge 
field according to 

s^^l d^xyf^l - 7^ + 2^a^<^!)''a^(/)" 
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where ga ^ 0, since the problem is not significantly more 
complicated. Rather than using the coordinates which 
were exploited in |l5j , we will take static spherically sym- 
metric solutions with coordinates 
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where U — U{r) and R = R{r), and dQ2 is the standard 
round metric on a 2-sphere. We take R > without 
loss of generality. We are primarily interested in black 
holes solutions with at least one regular horizon which 
are asymptotically flat as R ~f oo. 

If we assume that the gauge field is given by a magnetic 
monopole configuration with components 
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in an orthonormal frame then the Maxwell-type equa- 
tion. 
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is satisfied identically, and the remaining field equations 
obtained from variation of the action Q take the form 




Qa = ^aQ 



(10) 



and a prime denotes d/dr. By virtue of the Bianchi iden- 
tity one of Eqs. inj)-® can be derived from the others. 

One should note that in contrast to many simpler mod- 
els 0,0,^3 there is no simple duality relation between 
magnetic and electric solutions in the theory, on account 
of the multi-scalar exponential coupling term which mul- 
tiplies the electromagnetic part of the action Q. Thus 
we cannot simply read off the properties of solutions with 
an electric field in place of the monopole ansatz Q). It 
might be tempting to think of the action as a special 
case of an alternative model in which each scalar field is 
associated with an independent U{1) gauge field. How- 
ever, Q is not a simple truncation of such a model, since 
such an increase of the number of U{1) fields would lead 
to A'^ independent Maxwell-type equations, each with its 
own single scalar coupling, rather than a single equation 
of the form (O . Thus the present model can be expected 
to display differences in comparison to models in which 
moduli and several independent U{1) fields are present, 
allowing for the choice of duality-preserving combina- 
tions liilii. 



A. Nature of horizons 

A straightforward proof by contradiction of the type 
used in Refs. may be used to establish that the solu- 
tions of possess at most one regular horizon. We 
first note that the sum of Eqs. iQ and © gives 



which may be integrated to yield 
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for arbitrary constants a and (3. We assume that (|12f) 
possesses at least one real zero in order that there may 
exist at least one horizon and thus rewrite (|12|l as 



{r -r-){r -r+), 



(13) 
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FIG. 1; The scalars $ and ^, solutions of the model (1) with 
X2 — 1, q — 1, plotted for fixed values of r+ and r_, and 
variable third independent parameter. Case (a) corresponds 
to T,<s> = E>s, and the other cases are ordered according to 
decreasing ratio E4, /E* . 



where a = r+ + r- and f3 = r^r^, and we may assume 
without loss of generality that r_ < r-|-. We substitute 
(221 into to obtain 



[(r - r_)(r - r+)(?i" 
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Let us now suppose that, r_|_ 7^ r_, both values r = r± 
correspond to regular horizons, and that (?a > for any 
one of the scalar fields. Since the scalar is smooth at the 
horizon, if we evaluate (|14|l at r = r+ we obtain 
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from which it follows that (l)°''{r+) < 0. Similarly, if we 
evaluate (|14() at r = r_ we see that 0°'(r_) > 0. Now 
given that (r) is assumed to be smooth, it follows that 
it must have a maximum at an intermediate value r = vq 
such that r_ < ro < r+. However, if we evaluate ifTljl at 
r = To we obtain 



iro~r_)iro~r+)r"iro) 
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(16) 

from which it follows that (/)° (ro) > 0, giving a mini- 
mum, which is a contradiction. Thus if (p'^ is regular at 
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r = r_|_ it cannot also be regular at r = r_. The point 
r — r_ should thus correspond to a curvature singularity. 

Similarly, if we assume that ga < 0, then all the signs 
in the above arguments are reversed but we still obtain 
contradiction. It therefore follows that the solutions can 
at most possess one regular horizon, at r = r+. 

It is also useful to note that if ga > then 0" must be 
monotonically decreasing on the domain of outer com- 
munications of regular black hole solutions, since if it 
reached a minimum at a finite value rg > then at 
such a point Eq. IjlSfl would once again be true, but now 
with the implication that (p°'"{ro) < 0, again giving a con- 
tradiction. Likewise if < then 0" is monotonically 
increasing for r > r+. 

The function R(r) must be monotonically increasing 
in the domain of outer communications for either sign 
of ga, since by ijHl R" < at finite r for any solutions 
with non-trivial scalars. This leaves a global maximum 
as the only possible turning point for the function R{r), 
but such a choice would be inconsistent with asymptotic 
flatness, since by (|13|l R^r as r ^ oo, given e^^ 1. 



B. Constraints on charges 

An additional first integral of the field equations may 
be extracted as follows: take the difference of (0 and ((HJ, 
eliminate terms involving (e^^) using (|13|l . and elimi- 



nate terms involving e" 
resulting equation one obtains 



using ©. Integrating the 



N 



^ 9a 



'e^^ -r + c 



(17) 



a=l 



where c is an arbitrary constant. In order to obtain solu- 
tions which are regular at the outer horizon, it is in fact 
necessary to choose c — r+. With this choice, and again 
using ifT^ . ((T7jl may be integrated to yield 
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where is an arbitrary constant. From H18|l we see that 
r = T- will correspond to a singularity, as expected from 
above. 

Unfortunately, it does not appear to be possible to ob- 
tain an analytic solution in closed form to the remaining 
field equations. Since one of the N scalar equations © 
can be eliminated with the use of IjlHII . and since the 
function e^^ is given in terms of the function R accord- 
ing to H13|l . we are left with one first order ODE for R, 
namely 

(— + (r~r+)(r-r_)-2-(r-Af) + l 
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coupled nonlinearly to iV — 1 independent second order 
ODEs 114|) for the scalars. This is equivalent therefore 
to a system of 2N — 1 first order ODEs, which can be 
solved numerically. 

Much useful analytic information about the solutions 
can be obtained in relation to the values of the ADM 
mass, M, and the N scalar charges, Sa, which corre- 
spond to the 0(r~"'^) terms in the asymptotic expansions 
at spatial infinity. 
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The 0(r) coefficient, of the function R^ is a gauge 
parameter whose value fixes the choice of origin of the 
radial coordinate, r. We will use this gauge freedom to 
set 

i?i = 0. 

Expanding Eq. H18|l at spatial infinity by use of l|21l) and 
(|22(l it follows from the leading order term that the con- 
stant is related to the moduli vacuum charges, 0^, 
according to 
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Furthermore, if we also make use of H2U|I it follows from 
the next to leading order terms in H13|) and H18|) that the 
following relations hold between the constants r± and the 
asymptotic charges 



r± = M± \ M 




(24) 



The constraint that r+ > r_ then yields the inequality 
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which is saturated for the extreme solutions for which the 
horizon is degenerate with the inner singularity. 

With definitions of the asymptotic charges in hand we 
can now integrate various field equations between the 
horizon, r = r+ , and spatial infinity to obtain constraints 
on the charges. If the scalar equations © are integrated 
on this interval, for example, we find that 



9aQa^ 
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We note that for solutions which are regular at the hori- 
zon, given that the integrand of (|26|l is positive it follows 
that 



> 



a=l 



9a 



(27) 



5 



for each scalar charge, and furthermore Ea = if and 
only if Q = 0. Thus the charges E^ do not constitute an 
elementary primary scalar hair according to the defini- 
tion adopted in the Introduction. 

At first sight one might be tempted to assume that 
Eqs. H26(l provide constraints on all N scalar charges, 
and that we are therefore dealing with a system with 
purely secondary scalar hair. However, this is not in fact 
the case since the lower limit of integration, r+, already 
depends on the ADM mass and scalar charges according 
to if^ . 

The nature of the relations H2()|l is more readily un- 
derstood if we rewrite them in terms of functions 0" = 
0" — (f)'^, a = 1 . . . N, which have the leading order be- 
haviour (j)°'(r) ^ 'Sa/r + . . . at spatial infinity. We then 
see that equations (I26II can be rewritten as 



1 
'2ga 



In 
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for Eq 7^ 0. Since the bounds of integration are inde- 
pendent of the moduli vacuum charges, we see that the 
relations H2t)|l or (|28|l are constraints which determine the 
(j)'^ in terms of Qj and Eq. 

Let us turn to the question of whether there exist any 
constraints on the scalar charges Eq. In fact, there ap- 
pears to be only one additional constraint on the asymp- 
totic charges. This may be determined by noting that 
if one multiplies each of the scalar equations 10 by 
2B?'e^^ (jf'' and then takes the sum of the resulting N 
equations plus ^R^ (s^*^) times the difference of Eqs. 
© and 1)11(1 . one obtains the expression 
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We may integrate Eq. H29|) from r — r+ to spatial infinity, 
and use H13|) and (|24|l to obtain the following constraint 
on the charges 
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where 
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with 



Qa = Aae-S""^" 



(30) 



(31) 



(32) 



moduli vacuum charges. On account of the possibil- 
ity of different vacuum moduli charges, the individual 
scalars can effectively "see" different magnetic monopole 
charges, Qa- 

Using Eq. l(Hn|) we obtain the following expression 
equivalent to ((24|1 
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and is bounded above according to 
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The constraint (|3U|I reduces the number of indepen- 
dent scalar charges to — 1. Do any further constraints 
remain to be found? In the N = 2 case of two scalar 
fields this cannot be the case, given the numerical re- 
sults of Fig. 1 and the results of the dynamical systems 
analysis of ref. 15]: any further constraints would mean 
we no longer had a primary hair in contradiction with 
the results derived there. We will argue that no further 
constraints exist for N > 2 either. In particular, if fur- 
ther constraints exist on the charges then it would be 
possible to extract them from the field equations. If we 
consider the field equations at spatial infinity, then solv- 
ing order by order in inverse powers of r no constraints 
on the charges Ea are found, though we do find that all 
coefficients of terms 0(r~"), n > 2, in the asymptotic 
series (|2U|I - (|22() are completely determined. At the next 
order, for example. 
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Solutions with the asymptotic expansions (|20|I -H22 () in- 
clude many which correspond to naked singularities 
rather than black holes. The requirement that solu- 
tions also have a regular horizon leads to the further 
constraints if^ . l|Sn|) found upon integrating the inde- 
pendent field equations from r = r+ to spatial infinity, 
as above. However, we can obtain no more than one 
constraint for each independent field equation (|14|l . H19|l . 
and the relations (|28|l . (|30|l which give one constraint for 
each equation, exhaust the possibilities. Thus we find 
that there are A^ — 1 independent parameters among the 
A^ scalar charges, Eq. 



C. Thermodynamic quantities 



The quantity Q can be thought of as the magnetic 
monopole charge normalized by the weighted sum of the 



Even in the absence of complete analytic solutions, 
some thermodynamic relations may be obtained, given 
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that the black hole temperature and entropy are de- 
fined at the horizon, r_|_, which is related to the mass 
by 124|l . In particular, let us evaluate the derivative of 
(|13|l at the horizon. In terms of the surface gravity. 



i (e^^)' 



the horizon area A. 



r—r^ 



H 



47ri?2(r+), 



and using substituting for r_|_ from (|24() we then find 



M 



kA. 



H 



An 



N 



(38) 



We now define a magnetostatic potential, x(^)j according 
to 



drX = 
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We integrate (|40l) from r = 7-+ to r = oo, and use H2()|) 
to find that the magnetostatic potential of the horizon is 
given by 



1 ^ E 

^ a=l 



(41) 



It then follows that (|38|l is equivalent to the Smarr-type 
relation 
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III. LINEAR STABILITY 

In many cases hairy black holes are unstable. In this 
section we examine whether this is the case also for our 
solutions. It turns out that our solutions are stable 
against linear radial perturbations. The general pertur- 
bations of the solutions in the case of a single scalar field 
(dilaton) has been studied in great detail in |23|, using 
the methods of ref. and there stability has been 

proved. The calculations were, however, already very in- 
volved in that relatively simple case, and hence, following 
most of the literature on the subject [2ll l2^. we prefer 
to limit ourselves to the study of radial perturbations. 
Moreover, we shall consider only the case of two scalar 
fields. 

We consider the action (Q, and for simplicity we put 
q ^ 1, A2 = 1, since considering the more general set of 
coupling parameters in the N ^ 2 case of eq. ^ does 
not affect our conclusions. 

For the discussion of stability, it is convenient to use 
coordinates in which the metric takes the form 



= -e^(^-'^dt^ + eA(«'*)di?2 + R^dnl, (45) 



with 



<I> = $(t,i?), * = 1'(t,i?), (46) 



and the magnetic field is given in an orthonormal basis 
by 



PoJ.-Q ^0,0. (47) 
In these coordinates, the field equations read 



For completeness, we also note that according to (I24II . 
the mass M is a homogeneous function of degree ^ in 
Ay^, and of degree one in each Eq. The appropriate first 
law of black hole mechanics for the system is therefore 



N 



(43) 



where we have identified the temperature, T = K/{2n), 
and entropy, S — jA^, in the usual fashion. There is 
no independent variation dQ in (|43(l . since Q is related 
to the scalar charges via the constraint piljl . Indeed, for 
adiabatic variations 



dE„ 



(44) 



Interestingly enough, given the form of eq. H24|l for r+, 
there is no further contribution to the first law from in- 
dependent variations of the vacuum moduli charges, 
as there is in the case of other theories such as those 
of refs. [2^. In view of the relations (|28|l this is to be 
expected. 
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where the prime and the dot denote differentiation with 
respect to R and t, respectively. 

We perturb the field equations by time-dependent lin- 
ear perturbations of the form 

T{R,t) = T{R) + Sr{R,t)e^'^\ 

A{R,t) = A{R) + 6A{R,t)e''^\ 

<i>{R,t) = <i>{R) + S<l>{R,t)e^'^\ 

= + (5«'(i?,t)e*"*, 

where the perturbations are assumed small and the func- 
tions T{R), A{R), and denote the time- 
independent unperturbed solutions of the field equations. 
We did not perturb the Maxwell field since the electro- 
magnetic Bianchi identities imply that the monopole-like 
solution H47(l must be independent of the radial coordi- 
nate. 

The perturbed equations then read 
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boundary conditions it yields 

5A = 2i?($'<5$ + *'(5^'). 



(60) 



The problem of stability can then be reduced to the 
study of the perturbation of the scalar fields $ and ^ pH 
[2^ . After long manipulations of the perturbed equations, 
one can obtain a coupled system of second order linear 
equations for 5$ and (5^: 
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r- A' 
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where 
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-2i? 



r'-A' 



+2 



2 i? 

r'-A' 2 , , 
*' 



/2 



C{R) 



-2R 



2 R 
$'15" + $"v|/' 



(64) 



r'-A' 3 , , , 



(65) 



In Schwarzschild coordinates (|45ll the previous equa- 
tions are not regular at the horizon. Therefore, it neces- 
sary to define new "tortoise" coordinates, given by 

R*= I e(r-A)/2dE. 



Defining new fields 

M = i?(5$, v = R5^, 

and using the explicit time-dependence of the perturba- 
tive modes, one can finally put the stability equations 
(|61|l - (|62f) in the Schrodinger form 



di?*2 



?\i = Au, 



(66) 



where u is the vector of components (u, v) and A is a 
symmetric matrix with entries 



= -4($'(5$' + *'(5*') 


+ ^[c^-2*(5A-25$) 


^11 = -{a 


r'-A' 

2R 


+c^-2*(5$_2(5^') 


(59) 


A12 = A21 = 








A22 = -{b 


r'-A' 


Eq. H58|) can be immediately integrated. With suitable 


2R 



,r-A 



,r-A 
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The matrix A can be diagonalized, and has eigenvalues 



The classical solutions are stable under linear perturba- 
tions if the potentials Vi,2 are everywhere positive. This 
can be proved by generalizing the arguments of Chan- 
drasekhar |2^. In fact, Ht)6|) can be written as 



92 u 



Au = 0. 



(67) 



Multiplying I|fi7(l by dw^ jdt and integrating over i?*, one 
gets 



dt dt^ dt dR 



*2 



dt 



Au dR*. 



After integrating the second term by parts, and adding 
the complex conjugate equation, one obtains the energy 
integral 





du 


2 


du 


2 \ 


/( 


+ 


+ ut Au 




dt 




dR* 



const. 



(68) 



If the last two terms in the integral are bounded and 
positive definite, it follows that the integral of \du/dt\'^ 
is also bounded, ruling out any exponential growth of 
the perturbations. In our case, it is sufficient to show 
that u^Au is positive. This can be easily checked by 
diagonalizing the matrix A. If the eigenvalues Vi^2 are 
non-negative functions, then u^^Au is clearly positive. 

This can be readily checked for the exact solutions. In 
the other cases, one has of course to resort to numerical 
calculations. Let us consider for example the solution 
with $ = ^ [la. This can be written as 



d.^ = Jjl^Illi^^I^dt' 



,4/3 



j.i/3 (r-r+)(r-r_)4/3 
+r2/='(r-r_)4/3dl72, (69) 



(-t) 



2/3 



(70) 



In terms of the coordinate R such that the metric takes 
the form 14511 . one has 

_ (a2 - 2qA + A2)V3(a - bA -|- A2) ^^^^ 



(a2 -h A-h A2)4/3 



i?(a2- A2)2(a2+aA + A2)i/3 



(i?3 + 4a3)A2(a2 - 2aA + A^y/^a^ - bA + A^) ' 

(72) 



-2*_ -2* /a2-2aA + A2^'/3 



a'^ + aA + A2 



(73) 



0.2 



0.15 




0.05 
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FIG. 2; The potentials Vi and V2 for the exact solution, with 
different values ofb/a, and hence of the horizon. The extremal 
case corresponds to horizon at R — 0. 



where a — r_/3, b - 



r_/3, and 



i?3 + 2a3 + (74) 

In these coordinates, the singularity is located at i? = 0, 
and the horizon a,t R = {b + af/'^{b - 2a)^/^. 

One can now substitute the metric functions in H66|l . 
The analysis is greatly simplified by the fact that for this 
solution <& = ^1 and hence A{R) — B{R). The equations 
readily separate into two independent equations for u -|- ?; 
and u — V, with potentials Vi = An + A12 and V2 = 
All — A12, respectively. 

The potentials are plotted for different values of the 
ratio r^/r- in fig. 2. They vanish at the horizon and 
at infinity and are regular and positive in the interval. 
We can hence deduce the stability of the solution with 
$ = 

In the general case, numerical calculations show that 
the behaviour of Vi^2 is qualitatively the same as for exact 
solutions. Some examples are given in fig. 3. We can 
conclude that all the classical solutions are stable against 
radial linear perturbations. 



IV. DISCUSSION 

We have generalized the solutions previously found in 
ref. ^3 for magnetically charged black holes holes non- 
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Vi 




20 40 60 80 100 120 



V2 




R 



20 40 60 80 100 



FIG. 3: The potentials Vi and V2 for numerical solutions with 
fixed values of r+ and r_ , and variable third independent pa- 
rameter. The height of the peak increases with the value of 
the ratio between the scalar charges. The lowest one corre- 
sponds to E4, — E*. 

minimally coupled to two scalar fields, to the case of 
multiple scalar fields non-minimally coupled to a single 
magnetic monopole. Even though the complete analytic 
solutions have not been derived, we have succeeded in 
integrating enough of the field equations that constraints 
on the masses and charges can be derived. This analysis 
supports the claim made in ref. 15] that the solutions 
possess a primary hair. In the case of N scalar fields 
there are N —I independent parameters among the scalar 
charges. 

We have further shown that in the case of two scalar 
fields, that the black hole solutions are classically stable 
to radial perturbations, a feature which is absent in the 



case of a number of other hairy black hole solutions. We 
have no reason to expect that the case of multiple scalar 
fields as given in Sec. II would be any different. 

Our present analysis indicates that the primary hair 
of the multiscalar black holes has quite novel features 
as compared to the case of other hairy black holes. In 
particular, while the black hole is characterized by new 
independent degrees of freedom which are defined at spa- 
tial infinity, these charges must necessarily vanish if the 
magnetic field is turned off. To distinguish such scalar 
hair from the case of elementary primary scalar hair, 
which would theoretical exist even in the absence of gauge 
charges, we identify this new form of hair as contingent 
primary scalar hair. This suggests a further refinement 
of the no hair theorems by the statement: In theories sat- 
isfying the DEC, black holes do not possess elementary 
primary scalar hair. 

Our analysis of the thermodynamic properties of the 
solutions has been limited to the derivation of a Smarr 
formula, and of the first law. A further step would be the 
thorough study of the thermodynamical properties of the 
solutions using explicit numerical solutions. 

We also remark that solutions with properties analo- 
gous to those of the model studied in this paper have 
been obtained in the case of Gauss-Bonnet gravity non- 
minimally coupled to two scalar fields [23|. It would be 
interesting to investigate if our results can be extended 
also to that case. 
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